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Abstract 

The close approach to the Earth of comet C/1996 
B2 (Hyakutake) in March 1996 allowed searches for 
minor volatile species outgassing from the nucleus. 
We report the detection of deuterated water (HDO) 
through its loi-Ooo rotational transition at 464.925 
GHz using the Caltech Submillimeter Observatory. 
We also present negative results of a sensitive search 
for the J(5-4) line of deuterated hydrogen cyanide 
(DCN) at 362.046 GHz. 

Simultaneous observations of two rotational lines of 
methanol together with HDO in the same spectrum 
allow us to determine the average gas temperature 
within the telescope beam to be 69 ± 10 K. We are 
thus able to constrain the excitation conditions in the 
inner coma and determine reliably the HDO produc- 
tion rate as (1.20 ± 0.28) x 10 26 s” 1 on March 23-24, 
1996. Available IR, UV and radio measurements lead 
to a water production rate of (2.1 ± 0.5) x 10 29 s -1 
at the time of our HDO observations. The resulting 
D/H ratio in cometary water is thus (29 ± 10) x 10~ 5 , 
in good agreement with the values of (SO.Sl^;®) x 10” 5 
(Balsiger et ai 1995) and (31.6 ± 3.4) x 10 -5 (Eber- 
hardt et ai 1995) determined in comet P/Halley from 
in situ ion mass spectra. The inferred 3 a upper limit 
for the D/H ratio in HCN is 1%. 

Deuterium abundance is a key parameter for study- 
ing the origin and the early evolution of the Solar Sys- 
tem and of its individual bodies. Our HDO measure- 
ment confirms that, in cometary water, deuterium is 
enriched by a factor of at least 10 relative to the pro- 
tosolar ratio, namely the D/H ratio in H 2 in the prim- 
itive Solar Nebula which formed from the collapse of 
the protosolar cloud. This indicates that cometary 
water has preserved a major part of the high D/H 
ratio acquired in this protosolar cloud through ion- 
molecule isotopic exchanges or grain-surface reactions 
and was not re-equilibrated with H 2 in the Solar Neb- 
ula. However, there are strong presumptions that in- 
terstellar grains were partly mixed in the early nebula 
prior to cometary formation with water reprocessed in 
the warm inner part of the nebula and transported by 
turbulent diffusion. Scenarios of formation of comets 
consistent with these results are discussed. 

Subject headings: comets, deuterium, Solar Nebula, sub- 
millimeter spectroscopy 


1 . Introduction 

Deuterium is a sensitive indicator of the origin of 
water ice in comets. Comet ary ices may have con- 
densed in the primitive Solar Nebula or may have 
come directly from the protosolar cloud from w T hich 
the nebula was formed. In both cases, deuterium frac- 
tionation must have occurred at low temperatures be- 
tween molecular hydrogen (the main reservoir of deu- 
terium) and water, so the D/H ratio in H 2 0 should be 
higher than- the D/H ratio in H 2 in the Solar Nebula 
(usually reported as the protosolar ratio), estimated 
to be ~ 3 x 1CT 5 (see Section 4.2). 

The difference between the two scenarios comes 
from the kinetics of isotopic reactions: in cold dense 
interstellar clouds, the fractionation occurs through 
ion-molecule reactions and grain-surface chemistry 
and is very fast, so that the equilibrium can be 
reached even at a low temperature, resulting in high 
D/H enrichment factors in deuterated species (Wat- 
son 1974; Irvine and Knacke 1989). By contrast, the 
deuterium fractionation in the Solar Nebula occurs 
through reactions between neutral species (e.g. Geiss 
and Reeves 1981). In this case, the rate of isotopic 
exchange depends on temperature, pressure and time. 
It goes to zero at the low temperature and pressure 
values expected in the outer part of the nebula. In 
addition, the fractionation stops when the hydrogen 
escapes from the nebula. Therefore, a deuterium en- 
richment higher than a factor of 3 with respect to the 
protosolar value is not expected in the Solar Nebula. 

In summary, a high (much higher than 3) D/H en- 
richment in cometary water would thus be consistent 
with an interstellar origin of the ices from which these 
objects formed. On the other hand, a low D/H ra- 
tio (close to the protosolar ratio) would be consistent 
with a turbulent and well mixed nebula, at least out 
to the region of cometary formation. 

Until recently, the D/H ratio measured in the wa- 
ter of P/Halley was too uncertain to permit firm con- 
clusions (Eberhardt et ai 1987). However, the re- 
sults from the mass spectrometers aboard Giotto have 
been improved. The detailed analysis of the data 
of the ion mass spectrometer (IMS) by Balsiger et 
ai (1995) and those of the neutral mass spectrome- 
ter (NMS) by Eberhardt et ai (1995) independently 
lead the two groups to precise and consistent esti- 
mates of D/H in water of (SO.St^) X 10~ 5 (IMS) 
and (31.6 ± 3.4) x 10~ 5 (NMS). Such a high value 
implies an interstellar origin for at least a substan- 
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tial part of the P/Halley material. However, it might 
be argued that this comet may be atypical and not 
representative of most comets. It is thus valuable to 
measure the D/H ratio in other comets. 

The apparition of the bright comet C/1996 B2 
Hyakutake provided an opportunity to detect HDO 
from the ground by means of heterodyne spectroscopy 
at submillimeter wavelengths. Indeed, soon after 
its discovery on January 30, 1996 (Nakamura and 
Nakano 1996), this long-period comet turned out to 
be an intrinsically bright comet destined to make a 
very close approach to Earth on March 25, 1996 at 
only 0.102 AU (Marsden 1996). Due to the excel- 
lent viewing geometry, this comet proved to be an 
ideal and unprecedented target. Observations of the 
HDO loi-Ooo line at 464.925 GHz were organized on 
short notice at the Caltech Submillimeter Observa- 
tory (CSO). Good atmospheric conditions allowed an 
unambiguous detection (Lis et al 1996a) . This is the 
first detection of cometary HDO from the ground. In 
addition a sensitive search for DCN w r as performed. 

Section 2 describes the observations. They are 
part of the overall CSO observational campaign for 
comet C/1996 B2 (Hyakutake) presented by Lis et al 
(1997). In Section 3, the HDO production rate and 
HDO/H 2 O abundance ratio are derived. Section 4 
discusses the D/H ratio and its implications for the 
origin of comets and the structure of the primitive 
Solar Nebula. The excitation model for HDO is pre- 
sented in Appendix A. 

2. Observations 

Observations of HDO and DCN in comet C/1996 
B2 (Hyakutake) were carried out using the 10.4-m 
Leighton telescope of the Caltech Submillimeter Ob- 
servatory located on Mauna-Kea, Haw*aii. A log of 
the observations is given in Table 1. Other CSO ob- 
servations of comet Hyakutake are presented by Lis 
et al (1997). 

The loi -Ooo line of HDO at 464.92452 GHz (DeLu- 
cia et al 1971) was observed on March 23.5 and 24.5 
1996 UT, when the comet was at 0.11-0.12 AU from 
the Earth and 1.1 AU from the Sun. The observa- 
tions were performed near meridian transit when the 
comet was at 60-80° elevation. Observing conditions 
were quite good, characterized by zenith opacities at 
225 GHz near 0.05-0.06, as measured by a tipping 
radiometer. The single sideband (SSB) system tem- 
perature referred to above the Earth’s atmosphere at 


the elevation of the comet was about 1500 K at 465 
GHz. 

Observations of the J(5-4) transition of DCN at 
362.04648 GHz (Pickett et al 1996) were carried out 
on March 26.5, 1996 UT. The data were taken under 
good weather conditions characterized by a 225 GHz 
zenith optical depth of 0.06-0.07. However, the SSB 
system temperature was relatively high (—1200 K) 
due to the low elevation of the comet (27-30°). 

Comet Hyakutake was tracked using an ephemeris 
kindly supplied by Don Yeomans (private communi- 
cation). The ephemeris used on March 23 and 24 is 
based upon orbital elements referred to as Yeomans’ 
“Solution 18” which were computed from 354 astro- 
metric positions obtained from Jan. 1 to Mar. 21, 
1996. “Solution 21” was used during the DCN obser- 
vation on March 26. Ephemerides and pointing were 
checked frequently by performing five-point maps of 
HCN transitions and found to be accurate to ~3" 
r.m.s. 

The data of March 24 and 26 were acquired in 
the beam-switching mode, by chopping the secondary 
mirror at a 0.8 Hz rate with a throw of 4' in azimuth. 
Those of March 23 were made in position-switching 
mode, with the reference position offset by 10 The 
observations were made using the facility SIS receivers 
operating in a double-sideband (DSB) mode with sig- 
nal and image bands separated by 3 GHz. In order to 
minimize the system temperature, the receivers were 
tuned to the HDO and DCN frequencies in the upper 
sidebands. We used simultaneously two 1024-channel 
acousto-optical spectrometers (AOS) as backends: a 
low-resolution AOS which covers a total bandwidth of 
560 MHz with an effective spectral resolution of —1 
MHz (0.6 km s _1 at 465 GHz) and a high-resolution 
AOS of 50 MHz bandwidth providing a spectral res- 
olution of -100 kHz (0.06 km s _I at 465 GHz). 

The data were calibrated using the standard cho- 
pper-wheel method to give antenna temperatures out- 
side the atmosphere. The intensity calibration w*as 
checked by observing IRC +10216. The calibration 
uncertainty is typically 10% for lines observed with 
the 50 MHz spectrometer. For the interpretation of 
the observed signals, antenna temperatures were con- 
verted into main beam brightness temperatures using 
the relation Tb = T^/rimb, with efficiencies r/ mh mea- 
sured on Jupiter (Table 2). The telescope beam size 
(HPBW) at the distance of the comet corresponded 
to approximately 1200 km for HDO observations and 
1600 km for DCN observations. 
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Fig. 1. — The loi-Ooo line of HDO at 464.925 GHz in C/1996 B2 (Hyakutake) observed on March 23.5-24.5 UT 
with the 50 MHz spectrometer. The velocity scale is with respect to the comet nucleus. The spectrum has been 
smoothed to a spectral resolution of 0.5 km s” 1 . 


The DCN line was not detected. HDO was de- 
tected independently on March 23.5 and 24.5 UT 
(Table 3). The signal-to-noise ratio is ~7, when all 
observations are averaged. Figure 1 shows the aver- 
age spectrum obtained with the high-resolution spec- 
trometer. The average HDO spectrum acquired with 
the low-resolution spectrometer is shown in Fig. 2. 
Two methanol lines are present in the spectrum: the 
(9,2)-(9,l) A+-A~ line at 464.83467 GHz and the 
(15,0)-(14,1) A+ line at 461.75360 GHz in the image 
band (frequencies from Xu and Lovas 1997). The ob- 
served integrated line intensities of HDO and DCN 
are given in Table 3. They were measured on the 
high-resolution spectra in the -1.5 to 1.5 km s^ 1 ve- 
locity interval. Lines of other molecules observed with 
a high signal-to-noise ratio during the same period 
were found to be confined to this velocity range (Lis 


et al 1997). 

3. D/H ratio in water and hydrogen cyanide 

3.1. HDO production rate 

The HDO column density TVhdo is related to 
the observed line intensity integrated over velocity 
/ T B dv by: 

f hc z 

J Tedv = s^ r l AulfuNmo (1) 

where A u i is the Einstein ^-coefficient of the transi- 
tion, v ul is the rest frequency, and f u is the fractional 
population of the upper level. Eq.(l) assumes that 
the emission is optically thin and that the excitation 
conditions do not vary significantly within the region 
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Fig. 2. — Spectrum of the loi~0oo line of HDO at 464.925 GHz in C/1996 B2 (Hyakutake) observed on March 
23.5-24.5 UT with the 560 MHz spectrometer. Two CH 3 OH lines are present in the spectrum: the (9,2)-(9,l) 
A + -A~ and the (15,0)-(14,1) A + lines at 464.8347 GHz and 461.7536 GHz, respectively. The spectral resolution 
is 1.08 MHz (0.7 km s -1 ). 


sampled by the beam. 

As shown in Appendix A and in Fig. 5, with a 
beam size of 1200 km in diameter, the observations 
are sensitive to the collisional region at thermal equi- 
librium (LTE), where the population of rotational lev- 
els is described by the Boltzmann distribution: 


/« = 


9u 


Z(T kin ) 




( 2 ) 


where g u and E u are the statistical weight and the 
energy of the upper level, respectively, T k in is the 
gas kinetic temperature, and Z{T k in ) is the partition 
function. In a first approximation, it is thus relatively 
straightforward to infer the HDO column density from 
the observed line intensity, as long as the gas kinetic 


temperature can be reliably determined. 

According to hydrodynamical models of cometary 
comae, the gas kinetic temperature is expected to 
vary as a function of the distance from the nucleus. 
The HDO molecules contributing most (> 80%) of 
the observed signal are situated within 100 to 2000 
km from the nucleus. The gas temperature given by 
theoretical calculations is not reliable enough, due to 
the number and complexity of heating and cooling 
processes, so direct temperature measurements have 
to be used. The two methanol lines present in the 
HDO spectrum (Fig. 2) allow a direct estimation of 
Tkin within the field of view\ Their relative intensities 
are consistent with a CH 3 OH rotational temperature 
of 69 ± 10 K (Lis et al 1997). Other temperature 
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measurements are also available (Lis et ai 1997; Biver 
1997) and are summarized in Table 4. They all agree 
with a rotational temperature of about 70 K. There 
is no significant variation of the rotational tempera- 
ture with beam size. Calculations with a full excita- 
tion model (Bockelee-Morvan et ai 1994; Biver 1997) 
show that CO and CH3OH are close to LTE in the 
sampled region, so that their rotational temperature 
follows the kinetic temperature of the gas. Adopting 
a kinetic temperature of 69 ± 10 K for interpreting 
the HDO observations is thus secure. In the LTE ap- 
proximation, an uncertainty of 10 K in the assumed 
temperature affects the population of the loi rota- 
tional level and the inferred HDO column density and 
production rate by -—16%. 

For a more precise evaluation, we have to take 
into account radiative decay and infrared pumping, 
which affect the average population of the loi rota- 
tional level. Using the excitation model presented in 
Appendix A, we found that the upper level popula- 
tion within the CSO beam is increased by 30-40% 
compared to that inferred from the LTE model at 69 
K, depending whether collisions with electrons are in- 
cluded in the computations or not. However, the pop- 
ulation of the loi level in the full excitation model is 
less sensitive to the temperature: a 10 K temperature 
variation affects the population by only ~11%. 

The HDO column density deduced from the com- 
prehensive excitation model is (9.2 ± 1.9) x 10 12 
cm -2 for March 23-24. In order to convert the col- 
umn density into a HDO production rate, we assumed 
a Haser density distribution. A more realistic angu- 
lar density distribution is not required. Numerical 
simulations show that production rates calculated in 
cases of strongly anisotropic patterns do not differ by 
more than 20% from those calculated in the isotropic 
approximation (Biver 1997; Biver et ai 1997). In 
the case of comet Hyakutake, the HDO line (Fig. 1) 
as well as lines of other molecules detected with high 
signal-to-noise ratios at CSO present only moderate 
velocity shifts (Lis et ai 1997). Maps of the HCN 
coma also suggest a roughly isotropic gaseous coma, 
although some deviations from spherical symmetry 
are present (Lis et ai 1996b). Thus, we do not ex- 
pect our HDO production rate calculations to be sig- 
nificantly affected by the isotropic approximation. 

Following Crovisier (1994), the photodissociation 
rate was taken to be 1.3 x 10“ 5 r^ 2 s -1 , where is 
the heliocentric distance. The HDO production rate 
is only weakly sensitive to this parameter. We as- 


sumed an expansion velocity of 0.75 km s -1 . This 
outflow velocity is consistent with the shapes of the 
radio lines observed in C/1996 B2 (Hyakutake) at the 
end of March 1996 (Lis et ai 1997; Biver 1997). The 
inferred HDO production rate scales proportionally to 
the assumed outflow velocity, so that a velocity 10% 
higher would result in a HDO production rate 10% 
higher. 

Finally, we tested the sensitivity of the derived pro- 
duction rate to pointing errors. The effect is rather 
small. An offset of 5" in the pointing (those effectively 
observed during pointing checks were 3-5" at most) 
would lead to underestimating the HDO production 
rate by 10%. 

In Table 5, we give the 1 a confidence limit in the 
HDO production rate, taking into account the sta- 
tistical uncertainty along with all the uncertainties 
discussed above. 

3.2. Deuterium abundance in water 

To determine with good confidence the D/H ratio 
in water (D/H)h 2 Oj the water production rate at the 
time of the HDO observations is needed. Table 5 sum- 
marizes determinations available from different tech- 
niques. Ironically, the H2O production rate seems not 
as well constrained by observations as the production 
rate of HDO. 

Water was directly observed on March 24.5 through 
several ro- vibrational lines in the vibra- 

tional band at 5084 cm -1 using the NASA Infrared 
Telescope Facility (IRTF; Mumma et ai 1996). A 
production rate Qh 2 o = (2.3 ± 0.3) x 10 29 s*" 1 was 
derived (Dello Russo et ai 1997). 

Other measurements of the water production rate 
come from observations of its photodissociation prod- 
ucts. The A 2 £+ - X 2 Ili (0,0) band of OH near 3090 
A observed with the International Ultraviolet Ex- 
plorer (IUE) gives <2h 2 o = 1.6 - 1.9 x 10 29 s -1 (Fes- 
tou et ai 1997) for the March 23.1-25.1 period. The 
same band monitored from the ground by narrow- 
band photometry (Schleicher 1997) yields Qh 2 o = 
2.5, 2.7 and 1.7 x 10 29 s” 1 on March 23.5, 24.5 and 
25.1 UT respectively, using model parameters per- 
taining to the OH collisions! coma sampled by the 
instrumental aperture (Table 5). OH was also ob- 
served in the radio from its maser lines at 1667 and 
1665 MHz (Gerard et ai 1998), from which Qh 2 o was 
estimated to be ^ 2 x 10 29 s" 1 on March 24.1, and 
^ 1.4 x 10 29 s -1 on March 26.1. From observation 
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of HI Lyman-a, emission with the SWAN instrument 
on board SOHO, Bertaux et al (1998) derived Qh 2 o 
~ 1.9 x 10 29 s -1 . Finally, Hicks and Fink (1997) 
reported ~ 3 x 10 29 s“* from observations of the for- 
bidden line of oxygen 01 ( l D) at 6300 A. However, 
this determination is hampered by the uncertainty in 
the branching ratio to Ol{ l D) following H 2 0 dissoci- 
ation: it will not be considered here. 

In summary, the water production rate, as esti- 
mated by various techniques, is between 1.6 and 2.7 
x 10 29 s” 1 during the period of our HDO observa- 
tions. Many of these results are of preliminary nature 
and further work may lead to a better agreement be- 
tween the different estimates. It is presently difficult 
to make a choice between the various measurements, 
because each value relies on state-of-the-art models 
and assumptions. Even the determination made from 
the direct observation of water by its v\ + z/ 2 + ^3 — v \ 
band should be taken with a grain of salt, because the 
strength of this band is unknown and the results rely 
on the assumption that it is identical to that of the 
V 2 + vz band. 

The various data sets show significant temporal 
variations in the water production rate during the 
March 20 to 25 period (reaching 50% or more), which 
may be linked to the disruption of the nucleus ob- 
served in the visible (Jorda et al. 1996). Most agree 
with a surge of outgassing associated with the frag- 
mentation, followed by a plateau or a decrease in 
activity (Bertaux et al 1998; Gerard et at. 1998; 
Schleicher 1997; Table 5). Such temporal variations 
were also seen in the production rates of other species, 
such as HCN (Lis et al 1997). In addition, a peri- 
odic modulation of ~ 20% amplitude, related to the 
6.2 hours nucleus rotation, is observed in the OH gas 
production curve monitored by narrow-band photom- 
etry (Schleicher 1997; Schleicher et al 1998). The 
calculation of the £?hdo/Qh 2 o ratio thus requires a 
comparison of the amounts of HDO and H 2 0 emitted 
at the same time. The comparison between various 
data sets is not easy due to the difference in the fields 
of view (from a few arcsec to several degrees) and the 
time lag between production of HDO and H 2 0 and of 
the water photodissociation products. We note that 
the decrease seen in the HDO line intensity between 
March 23.5 and 24.5 (Table 3), which is marginal tak- 
ing into account combined uncertainties, is neither ob- 
served in the OH data of Schleicher (1997) obtained 
with a similar field of view, nor in the CH3OH lines 
observed in the HDO spectrum (Lis al 1997). 


Using a water production rate of (2.1 ± 0.5) x 10 29 
s -1 (average of the values given above) together with 
the average HDO production rate for March 23-24 
of (1.20 ± 0.28) x 10 26 s" 1 , we infer a Qhdo/Qh 2 o 
ratio of (57 ± 19) x 1CT 5 , wdiich corresponds to a 
(D/H)h 2 o ratio of (29 ± 10) x KT 5 . 

3.3. Deuterium abundance in hydrogen cya- 
nide 

The 3 a upper limit for the DCN J(5-4) line inten- 
sity on March 26 (Table 1) corresponds to Qdcn < 
2.6 x 10 24 s" 1 . On the same date, the *7(4-3) line of 
HCN at 354.505 GHz was observed, from which Qhcn 
= 2.1 - 2.4 x 10 26 s” 1 was derived (Lis et al 1997). 
Both production rate determinations were made con- 
sistently, using the same model (Crovisier 1987) and 
model parameters. The upper limit for the (D/H)hcn 
ratio in Hyakutake is thus approximately 1%. 

4. Discussion 

The D/H ratio measured in water in comet Hyaku- 
take is consistent with the determination of D/H in 
water in comet P/Halley obtained by NMS (Eber- 
hardt et al 1995) and IMS (Balsiger et al 1995) 
aboard Giotto. It is also in agreement with the anal- 
ysis of the submillimeter detection of HDO in comet 
C/1995 01 (Hale-Bopp) which gives (33 ± 8) x 10" 5 
(Meier et al 1998a). Note that, on the basis of their 
orbital elements, all three of these comets are believed 
to originate from the Oort cloud. 

4.1. Comparison with the interstellar medium 

The D/H ratio in water in these comets appears to 
be enriched by a factor of at least 10 with respect to 
the protosolar value (Fig. 3 and 4). Geiss and Reeves 
(1981) pointed out long ago that high isotopic enrich- 
ments could not have occurred via reactions between 
neutral molecules in the primitive Solar Nebula; they 
showed that if a high enrichment of D/H were discov- 
ered in solid objects in the Solar System, it would im- 
ply that the enrichment took place in the interstellar 
medium (ISM) through ion-molecule reactions, prior 
to the collapse of the protosolar cloud. More recently, 
it has been shown that high enrichments could also re- 
sult from grain-surface processes in the ISM (Tielens 
1983; Brown and Millar 1989). 

Owen et al (1986) were the first to present D/H 
measurements demonstrating that highly enriched 
components are present in the outer Solar System. 
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Fig. 3. — D/H in protosolar cloud and ISM today. Only the high D-enrichment observed in meteorites is drawn. 
The D/H ratio in cometary H 2 O measured in P/Halley (Eberhardt et al 1995; Balsiger et al 1995) is the value 
given without an error bar; that measured in C/1996 B2 (Hyakutake) is the value given with an error bar. The 
vertical arrow is the upper limit on the DCN/HCN ratio in comet Hyakutake. The sources of the other data are 
given in Section 4. 


Although subsequent observations have led to changes 
in the derived values (Section 4.2), the basic point re- 
mains: if condensed matter has never been able to 
equilibrate with hydrogen gas through thermal reac- 
tions and mixing in the Solar Nebula, it will retain the 
values of D/H set in the ISM. These values should dif- 
fer from one species to another, reflecting the different 
binding energies of the deuterated molecules. Figure 
3 compares cometary D /H ratios with those observed 
in various components of the ISM. Note that values 
of interstellar D/H’s represent the isotopic ratios in 
specific regions at the present stage of galactic evolu- 
tion, whereas the cometary D/H pertains to that of 
the protosolar cloud 4.55 Gyr ago. 

The D/H ratio in hydrogen in the local ISM has 
been investigated through observations of the Ly a 
lines of hydrogen and deuterium atoms in absorption 
toward several stars (Linsky et al. 1995; Wood et al 
1996; Piskunov et al 1997). The average D/H ratio 


is ~ 1.6 x 10" 5 with the individual values ranging 
from 1.4 to 2.2 x 10" 5 (Piskunov et al 1997). This 
is about a factor of 2 lower than the protosolar value, 
which reflects the local ISM deuterium abundance at 
the time and place of Solar System formation. 

HDO has been observed in several galactic hot 
cores (Jacq et al 1990; Helmich et al 1996; Gen- 
sheimer et al 1996). Hot core regions are parts of gi- 
ant molecular clouds where high-mass stars are form- 
ing. In such warm regions, some species observed in 
the gas phase are believed to be produced by the sub- 
limation of icy grain mantles (e.g. H 2 O, CH 3 OH), so 
that the gas composition of these species potentially 
traces the ice composition. HDO is observed through 
its loi - Ooo submillimeter line as well as several mil- 
limeter lines. The abundance of H 2 O is measured 
through the 3i3 — 220 line of H^O at 203 GHz. Most 
of the observed values correspond to (D/H)h 2 o in the 
range 10-30 x 10" 5 , comparable to the observed D/H 
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ratio in comet ary water. 

DCN has also been observed in hot cores, as well 
as in cold galactic clouds with no or minimal star for- 
mation. The deuterium enrichment in HCN is found 
to be even higher than in H 2 O with D/H close to 5-8 
xlO” 3 in hot cores, but as high as 2 x 10" 2 in cooler 
clouds as TMCl or the cores in p Ophiuchi (Wootten 
1987; Mangum et al 1991; Helmich et al 1996). Al- 
though the hot core values are just below the upper 
limit we measured for the DCN/HCN ratio in comet 
Hyakutake, our limit lies below the ratio observed in 
the coldest interstellar medium clouds. From the re- 
cent detection of DCN in C/1995 01 (Hale-Bopp), 
cometary (D/H)hcn is (2.3 ±0.4) x 10“ 3 (Meier et al 
1998b), i.e. significantly below the values observed 
in hot cores and cold clouds, although well above the 
protosolar D/H value. 

The chemistry necessary to produce such high lev- 
els of deuteration in the interstellar medium requires 
very low temperatures. The most current explana- 
tion given for the high D-enrichments observed in 
the warm (> 100 K) environment of hot cores is 
that D-rich molecules formed and froze onto grains 
in an earlier low-temperature phase. Reactions with 
atomic deuterium on grain surfaces may also con- 
tribute (Tielens 1983). Another more recent sugges- 
tion is that water ice mantles could be created in in- 
terstellar clouds by the passage of a C-type shock, fol- 
lowed by rapid hydrogenation of abundant gas-phase 
oxygen, and subsequent depletion of water onto the 
cold dust grains (Bergin et al 1998). This model does 
offer an alternative method to create water-rich man- 
tles. However, the level of deuterium fractionation is 
difficult to assess as a shock may alter the HD0/H 2 0 
ratio, as will any post-shock evolution prior to deple- 
tion onto the grain surface. 

In the hypothesis that grains infalling from the 
protosolar cloud onto the discoid nebula were never 
subsequently reprocessed or mixed with other compo- 
nents of the Solar Nebula (see discussion Sect. 4.3), 
the D/H ratios observed in comets could have inter- 
esting implications. From models of deuterium chem- 
istry in dense interstellar clouds (Millar et al 1989), 
the lack of a highly deuterium-enriched component 
in cometary HCN would suggest that the outer parts 
of the protosolar cloud were warmer than ~ 30 K 
when the gases condensed onto grains. These models 
give (D/H)h 2 o ratios consistent with those observed 
in comets for temperatures of ^ 50 K. However, as 
discussed in Section 4.3, there are strong presump- 


tions that the D/H ratios observed in comets are not 
unaltered relics of the protosolar chemistry. 

4.2. Comparison with other Solar System ob- 
jects 

A comparison of the deuterium abundance in co- 
mets with that in other Solar System bodies may have 
important cosmogonical implications. The D/H ratio 
in objects relevant to our problem, namely in objects 
containing relics of enrichments which occurred in the 
nebula at the birth of the Solar System, are shown 
in Fig. 4, significantly updated from Lecluse et al 
(1996). We do not include Mars and Venus which ex- 
hibit evidence of a strong atmospheric fractionation 
during the history of the planets, nor Titan since it 
was formed in the Saturn subnebula where chemical 
conversions and reprocessing of deuterium in carbon 
compounds presumably occurred (Prinn and Fegley 
1981; Lunine and Tittemore 1993). New data shown 
in this figure concern the protosolar value, meteoritic 
data, comet Hyakutake, and preliminary results for 
Jupiter and Saturn from the Infrared Space Observa- 
tory (ISO). 

The protosolar D/H ratio has been estimated by 
several authors using the conversion of D to 3 He in 
the early Sun to trace the original D from measure- 
ments of the 3 He/ 4 He ratio in the solar wind (Geiss 
and Reeves 1972; Geiss 1993; Niemann et al 1996). 
Gautier and Morel (1997) have recalculated the pro- 
tosolar D/H ratio by the same approach, but using a 
solar evolutionary model with microscopic diffusion. 
They derive (D/H ) pro to = (3.01 ± 0.17) x 10“ 5 . Geiss 
and Gloecker (1997) estimate that the helium isotopic 
fractionation results in lowering the 3 He/ 4 He ratio at 
the surface with respect to the solar wind value by 
1-9%, which leads to a lower limit for (D/H) pr0 * 0 of 
approximately 2.6 x 10“ 5 . Accordingly, we conserva- 
tively assume (D/H) pro to = (3.0 to a) x 10~ 5 . 

The terrestrial D/H reference value is currently 
15.6 x 10 -5 , measured in ocean water (SMOW) (De 
Witt et al 1980). However, the D/H ratio in oceans 
is expected to vary with time. Therefore, we prefer 
to adopt as a reference the averaged terrestrial value 
recently inferred by Lecuyer et al (1997) as equal to 
(14.9 ± 0.3) x KT 5 . 

Recent analysis of meteoritic data reveals the pres- 
ence of two different D-enrichments in water in the 
Semarkona and Bishunpur meteorites (Deloule et al 
1997, 1998). These authors find evidence of a minor 
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Fig. 4 — D/H ratios in the Solar System. The Earth value (Lecuyer et al 1997) and the protosolar value (see 
Section 4.2) axe shown for comparison. Asterisks indicate ground based observations. HD* for Jupiter is from 
Smith et al (1989b). The vertical arrows labelled HD* for Uranus and Neptune indicate upper limits derived by 
Smith et al (1989a) from searches of HD lines from the ground. The values deduced from HD observations with 
ISO are from Encrenaz et al (1996) (Jupiter) and Griffin et al (1996) (Saturn); the value obtained by Galileo in 
Jupiter is from Niemann et al (1996). The D/H ratio in methane given for Saturn and Uranus are from Noll and 
Larson (1990) and from de Bergh et al (1986), respectively. Two values deduced from the measured CH3D/CH4 
ratio are shown for Neptune, that derived from de Bergh et al (1990) from near-IR observations (large error 
bar) and that of Orton et al (1992) (small error bar). The D/H ratio deduced for protoneptunian ices (see text) 
is drawn as a rectangular box. The measurements in the Semarkona and Bishunpur meteorites, which show f the 
presence of two different D/H enrichments, are from Deloule et al (1997, 1998). For comet P/Halley, are shown 
the values of Eberhardt et al (1995) and Balsiger et al (1995). The value given for C/1996 B2 (Hyakutake) is 
from this work. 

(15%) highly enriched component with D/H = (73 ± 

12) x 10“ 5 and a major (85%) moderately enriched 
component with D/H = (8.8 ± 1.1) x 10“ 5 , namely 
less than the Earth value of 14.9 x 10" 5 . 

As initially recognized by Hubbard and McFarlane 
(1980), the observed D/H ratio in giant planets re- 
sults from mixing of hydrogen originating from the 


primitive Solar Nebula (and thus having a protosolar 
D/H ratio) with ices that may have been enriched in 
deuterium (for instance in the ISM) with respect to 
the protosolar value. In Jupiter and Saturn, the ice 
enrichment is currently considered negligible. How- 
ever, Guillot et al (1997) estimate that the amount 
of ices in Jupiter could be as high as 32 Earth masses 
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compared to 318 Earth masses for the whole Jupiter 
mass. In such a case, a detectable enrichment can- 
not be ruled out. For instance a D/H enrichment in 
proto-ices by a factor of 10 would result in an enrich- 
ment in Jupiter by 20% with respect to the protosolar 
value (Gautier and Morel 1997). 

Preliminary results from ISO based on the detec- 
tion of HD rotational lines on Jupiter and Saturn (En- 
crenaz et al 1996; Griffin et al 1996) are consistent 
with the protosolar value given above (Fig. 4). The 
preliminary results obtained by the mass spectrome- 
ter aboard the Galileo probe into Jupiter (Niemann 
et al 1996) give a lower limit also consistent with 
the protosolar value. The accuracy of this measure- 
ment is expected to be improved when the spectrom- 
eter will have been fully calibrated. At this point, 
we do not know whether the deuterium abundance 
in Jupiter is representative of the protosolar value or 
not. Note that the various estimates of the D/H ratio 
in Jupiter based on observations of CH4 are in conflict 
each other and spread out from 1 x 10 -5 (Gautier and 
Owen 1989) up to 5.5 x 10 -5 (Carlson et al. 1993; 
see Lecluse et al 1996 for a detailed discussion and 
other references). For clarity, these are not included 
in Fig. 4. 

The mass of hydrogen in Uranus and Neptune is 
much smaller than the mass of ices. If these planets 
were well mixed at least once in their history, as as- 
sumed in evolutionary models of Pollack and Boden- 
heimer (1989) or in the most recent models of their 
formation (Lissauer et al 1995), water and hydro- 
gen should be equilibrated at the high temperatures 
present in the interiors. In other words, the present 
hydrogen has acquired the deuterium enrichment ini- 
tially present in protoices embedded in the nebula and 
which formed the cores of Uranus and Neptune. The 
D/H ratio in protoices can be estimated (Lecluse et 
al 1996) from models of planetary interiors which 
give the relative proportions of ices, rocks and gas. 
Following Eq. (34) of Lecluse et al (1996), adopting 
the model of Podolak et al (1995) for the interior 
of Neptune, and taking the most precise determina- 
tion of D/H in the atmosphere, namely the value of 
Orton et al (1992) from CH 3 D, and (D/H ) pro to = 
( 3 . 012 : 5 ) X 10" 5 , we obtain a D/H ratio in protonep- 
tunian ices of (9.0lj;ij) x 1(T 5 . This value is com- 
parable with the moderately enriched component in 
the Semarkona and Bishunpur meteorites. Note that 
the recent measurements of the HD line at 37.7 pm 
in Neptune from ISO (Feuchtgruber et al 1997) is 


in good agreement with the Orton et al value. The 
present estimate of D/H in Uranus is uncertain, but 
an upper limit of 10 x 10” 5 derived from observations 
of HD lines in the visible suggests that proto-uranian 
ices were significantly less enriched in deuterium than 
cometary ices. 

In summary, two categories of fossil deuterium en- 
richment in water appear in the Solar System: (i) 
a moderate deuterium enrichment (a factor of ~ 3 
with respect to the protosolar value) is found in pro- 
toneptunian ices and the low enriched component of 
Semarkona and Bishunpur; (ii) a high deuterium en- 
richment is found in comets (a factor of ~ 10), and the 
highly enriched component of Semarkona and Bishun- 
pur (a factor of ~ 25). 

4.3. Implications for the origin of comets 

As discussed previously, the high deuterium en- 
richment in water can be obtained only through ion- 
molecule or grain-surface reactions occurring in the 
interstellar medium. However, an enrichment up to 
a factor of 3 can also be reached in the Solar Nebula 
(when neglecting diffusion) through thermal isotopic 
exchange (Lecluse and Robert 1994). Detailed cal- 
culations based on the generalization of the model of 
Dubrulle (1993) are out of the scope of this report and 
will be presented in a forthcoming paper (Drouart et 
al . 1998). It is assumed in this latter work that the 
transport of small icy particles and water vapor in the 
turbulent nebula is diffusive and that small particles 
remain mixed with the gas (Dubrulle et al 1995). It 
is also assumed that ices infalling from the protoso- 
lar cloud on the nebula disk are strongly enriched in 
deuterium, by a factor of about 25 as in the highly 
enriched component of the Semarkona and Bishun- 
pur meteorites. Preliminary results indicate that the 
diffusion occurring in the nebula that is required to 
transport the angular momentum by turbulent vis- 
cosity to Neptune redistributes the deuterium enrich- 
ment uniformly out to the orbit of this planet. Ac- 
cordingly, an enrichment of ^3 could have occurred 
at 3 AU where presumably meteorites were formed, 
as well as at 30 AU. But it is not possible to obtain 
an enrichment by a factor of 3 in one part of a tur- 
bulent nebula and by a factor of 10 in another part of 
the same turbulent nebula , except very briefly in the 
outer part of the early nebula when the mixing was 
not completely achieved . 

From these results, there are two possible scenar- 
ios, as detailed in Drouart et al (1998): 
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(i) Icy cometesimals might have been rapidly for- 
med — in less than 1 or 2 x 10 5 years — in the 
Uranus-Neptune distance range, as assumed in stan- 
dard models of cometary formation (Fernandez and 
Ip 1981; Duncan et al 1987). At this time, the mix- 
ing is not yet complete and the deuterium enrichment 
not yet uniform throughout the nebula. Comets of 
kilometer sizes are formed in a few hundred thou- 
sand years (Weidenschilling 1997). They are rapidly 
expelled towards the Oort cloud by early giant plan- 
ets. Uranus and Neptune continue to grow for sev- 
ered millions years (Pollack et al 1996), by accre- 
tion of planetesimals formed later than cometesimals 
from microscopic grains that axe mixed with the neb- 
ula gas and exhibit a deuterium enrichment of the 
order of 3, as in the low enriched component of the 
Semarkona and Bishunpur meteorites. This scenario 
requires an appropriate nebula, an appropriate tem- 
poral sequence for the formation and the expulsion of 
comets and a continuous formation of planetesimals 
from microscopic grains during the major part of the 
history of Uranus and Neptune. On the other hand, it 
permits us to explain the presence of crystalline sili- 
cates detected in several comets (Hanner et al 1994) 
and recently clearly confirmed by ISO observations 
of Hale-Bopp (Crovisier et al 1997). Hanner (1997) 
indicates a proportion of 20% of crystalline silicates 
with respect to amorphous silicates in comet Hale- 
Bopp. Silicates are observed in amorphous form in 
the ISM. Forming crystalline silicates requires tem- 
peratures higher than 1300 K, impossible to reach 
in the outer nebula. A partial diffusive transport of 
these species from the inner nebula out to the Uranus- 
Neptune region is a plausible solution of the problem. 

(ii) The nebula might have been turbulent out to 
the region of formation of Neptune, somewhere be- 
tween 25 and 30 AU, since several authors consider 
that this planet may have migrated outward by sev- 
eral AU during its formation process (Fernandez and 
Ip 1996; Malhotra 1998). Beyond Neptune, the turbu- 
lence may have vanished and a laminar regime became 
predominant, implying a low mixing with the turbu- 
lent part of the nebula. This scenario is consistent 
with observations of disks around young protostars 
which exhibit sizes generally larger (up to 1000 AU) 
than assumed in current models of turbulent nebula 
(Beckwith et al 1993; Strom et al 1993). Grains in- 
falling from the protosolar cloud on this outer nebula 
may have sublimated and condensed again, as dis- 
cussed by Lunine et al (1991), but the local temper- 


ature was low enough (Dubrulle 1993) that no iso- 
topic exchange with hydrogen could occur. Thus, 
even in this case, the ices preserved the deuterium 
enrichment acquired in the interstellar medium. In 
this scenario, comets P/Halley, Hyakutake and Hale- 
Bopp would have been formed in this part of the neb- 
ula, namely around the location of the present Kuiper 
Belt or beyond. Levison and Duncan (1997) pointed 
out that a substantial part of the objects originating 
from the Kuiper belt may have become part of the 
Oort cloud. Interpreting the presence of crystalline 
silicates in comets requires also some mixing of the 
laminar nebula with the turbulent one, as observed 
in laboratory experiments (Daviaud et al 1992) but 
not easily simulated numerically. 

Discriminating between these scenarios requires 
improvements in the theories of the origin of comets, 
as well as additional measurements of (D/H)h 2 o in 
both long-period and short-period comets. A low 
D/H ratio in Jupiter family comets originating from 
the Kuiper Belt would imply that the nebula was tur- 
bulent much farther from Neptune and would favor 
scenario (i). A D/H ratio of the order of or larger than 
that in P/Halley, Hyakutake and Hale-Bopp would 
favor scenario (ii) and a common origin of comets be- 
yond Neptune. 

Whatever the scenario of cometary formation, it is 
clear that cometary ices have an interstellar origin, as 
previously advocated by a number of authors (Green- 
berg 1982; Mumma et al 1993; Irvine et al 1996). 
However, both scenarios (i) and (ii) assume that in- 
terstellar grains were partly mixed during cometary 
formation with water reprocessed in the warm inner 
part of the turbulent nebula. This provides an in- 
terpretation of the presence of crystalline silicates in 
some comets and also explains why the (D/H)h 2 o val- 
ues observed in the three Oort cloud comets are some- 
what less than in the highly enriched component of 
the Semarkona and Bishunpur meteorites. 

4.4. Origin of water on the Earth 

The comparison of the D/H ratios in comets and 
in meteorites with that of the Earth provides clues 
for the origin of the oceans of our planet. The high 
cometary values make it difficult to defend a major 
cometary origin of terrestrial water. On the contrary, 
in whole rock carbonaceous meteorites and in chon- 
drules, the mean value of D/H, namely the average 
of the high enriched component and of the low en- 
riched component, is consistent, within uncertainties, 
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with the Earth value (Deloule et al 1997; Lecuyer 
et al 1997). The current assumption that plan- 
etesimals which formed the Earth were of chondritic 
composition favors the scenario of the formation of 
oceans mainly by outgassing from the Earth’s surface. 
Some water on the Earth could be inherited from a 
cometary bombardment subsequent to its formation 
(Owen and Bar-Nun 1995) but the deuterium abun- 
dances presently measured in comets suggest that this 
contribution should not exceed 10% (Deloule et al 
1997) to 30% (Owen 1997). 

5. Conclusion 

We have measured a (D/H)h 2 o ratio of (29 ± 10) 
x 10“ 5 in the water of C/1996 B2 (Hyakutake). We 
may expect in the future to reduce the uncertainty 
of this measurement, with better constraints on the 
water production rate of Hyakutake at the time of the 
HDO observations. 

This D/H ratio is consistent with the value found 
in P/Halley and more recently in Hale-Bopp (Meier 
et al. 1998a). The high enrichment found in these 
objects with respect to the protosolar value implies 
that ices which formed comets were initially enriched 
in the protosolar cloud through ion-molecule isotopic 
exchanges or grain-surface reactions, prior to the col- 
lapse of the cloud. 

However, we believe that some mixing with ices 
reprocessed in the turbulent part of the nebula subse- 
quently occurred due to diffusive processes. Indeed, 
the D/H enrichment found in these three comets is 
somewhat less than that found in some components 
of primitive meteorites. The presence of crystalline 
silicates in comets, which requires very high temper- 
atures achievable only in the inner part of the early 
nebula, favors this assumption. 

On the basis of the D/H enrichments in water 
found in other objects of the solar system, we dis- 
cuss two possible scenarios for the origin of comets 
which have strong consequences for the structure of 
the nebula. Following Drouart et al (1998), in both 
scenarios it is assumed that the transport of the an- 
gular momentum from the center of the nebula out to 
Neptune or farther was made by turbulent viscosity. 
Due to diffusion, the radial distribution of the deu- 
terium enrichment in water became rapidly uniform 
throughout the nebula out to Neptune at least. 

In the first scenario, comets P /Halley, Hyakutake 
and Hale-Bopp were formed very early in the his- 


tory of the nebula and rapidly expelled to the Oort 
cloud by growing giant planets, then keeping a part of 
the deuterium enrichment acquired in the protosolar 
cloud. Uranus and Neptune grew slowly during all the 
lifetime of the nebula by accreting planetesimals con- 
tinuously formed during the whole story of the nebula 
and exhibiting a low deuterium enrichment. In this 
scenario, comets coming from the Uranus-Neptune 
region populated the Oort cloud. 

In the second scenario, comets formed in the non- 
turbulent part of the nebula, in this outer region be- 
yond Neptune, which was only weakly mixed with the 
inner nebula. Icy grains kept most of their interstel- 
lar deuterium enrichment, whatever the time after the 
formation of the nebula. Comets were subsequently 
expelled by giant planets towards the Oort cloud. The 
deuterium enrichment in ices which formed Uranus 
and Neptune was low due to turbulent mixing. In 
this scenario, all comets in the Oort cloud originate 
from the Kuiper Belt or beyond. 

For the moment, we cannot discriminate between 
these two scenarios, since the D/H ratio in comets 
originating from the Kuiper Belt has not yet been 
measured. In addition, the three comets in which 
HDO has been observed may not be representative of 
the whole population of the Oort cloud. A low D/H 
ratio in Jupiter family comets originating from the 
Kuiper Belt would imply that the nebula was tur- 
bulent much farther from Neptune and would favor 
the first scenario. A D/H ratio of the order of or 
larger than that in P/Halley, Hyakutake and Hale- 
Bopp would validate the second scenario and favors 
a common origin of comets beyond Neptune. The 
presence of crystalline silicates in comets is consis- 
tent with both scenarios, provided that the mixing 
with the inner part of the nebula was large enough. 

We conclude that it is extremely important to mea- 
sure D/H in H 2 0 in every comet of sufficient bright- 
ness. Periodic comets may require visits by spacecraft 
and we look forward to the ROSETTA mission as the 
first such encounter. 
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APPENDIX 

A. Excitation and intensity of the HDO rota- 
tional lines 

Local thermodynamical equilibrium (LTE) is only 
valid in the inner coma where collisional excitation 
prevails. In the outer coma, the HDO rotational lev- 
els reach the fluorescence equilibrium resulting from 
the balance between solar excitation of vibrational 
bands and vibrational and rotational spontaneous de- 
cay. Modeling the full evolution of the excitation con- 
ditions of HDO in the coma is necessary to investigate 
to which extent observations sample the LTE region 
and how the rotational populations are affected by 
radiative decay and IR radiative pumping. Due to its 
strong dipole moment (/. ia = 0.657 D and /i# = 1.732 
D), HDO belongs to the category of molecules (which 
also includes H 2 O and HCN) whose rotational popu- 
lations relax rapidly toward fluorescence equilibrium 
where only the lower states are significantly popu- 
lated. 

The model computes the evolution of the rotational 
population distribution within the ground vibronic 
state as the HDO molecules expand in the coma, en- 
countering varying excitation conditions. This model 
is similar to those developed for e.g. H 2 0 (Bockelee- 
Morvan 1987), H 2 CO (Bockelee-Morvan and Cro- 
visier 1992) and linear molecules (Crovisier 1987). 
The rate equations include thermal excitation by col- 
lisions, absorption of the solar IR radiation in the 
vibrational bands and spontaneous decay. In addi- 
tion to collisions with water, the model also considers 
collisions with electrons (see Biver 1997, Biver et al 
1997). The main parameters of the model are the 
kinetic temperature of the gas, T k i n , the total cross- 
section for rotational de-excitation by collisions with 


H 2 0, <t c , and the outflow velocity v exp . The Haser 
model is used to relate the density of water in the 
coma to its production rate Q h 2 o* Details on the 
assumed electron density n e and temperature T e are 
given in Biver (1997) and Biver et al (1997). In sum- 
mary, the electron density as a function of distance 
to nucleus r follows the behavior of the total ion den- 
sity measured in P/Halley by ion mass spectrometers 
(Balsiger et al. 1986; Haberli et al 1995) with (i) are- 
combination zone in the inner coma at r < Res where 
the density varies according to (Qh 2 o\/ 77) 1//2 /r and 
the electrons are thermalized to the temperature of 
the gas; (ii) a pile-up region beginning at r = Res 
where the ions/electrons density increases as the re- 
sult of increasing electron temperature and decreasing 
electronic recombination (Eberhardt and Krankowsky 
1995); (iii) an external region at r > R rec where the 
density follows the Q/r 2 dependence expected for ra- 
dial expansion without significant additional ioniza- 
tion. Res and Rrec are scaled according to 

and respectively, and the absolute scaling is 

made such that the modeled n e reproduces the mea- 
surements in P/Halley for a water production rate of 
6 x 10 29 s -1 . With this model, Res = 3800 km and 
Rrec = 9000 km for comet Hyakutake at the time of 
the HDO observations. 

Energy levels, rotational line frequencies and stren- 
gths of the ground vibrational state were taken from 
the JPL molecular data base (Pickett et al 1996). 
For modeling radiative pumping of the vibrational 
bands, we used the spectroscopic data of the GEISA 
molecular data base (Husson et al 1992). In the 
present work we only consider the three fundamen- 
tal bands of vibration of HDO V\ (2727 cm -1 ), i/ 2 
(1402 cm” 1 ) and (3707 cm -1 ), whose total band 
strengths are 51, 209 and 129 cm -2 atm -1 , respec- 
tively. Harmonic and combination bands can be ne- 
glected on first order. This translates into solar ex- 
citation rates at 1 AU from the Sun of 6.5 x 10 -5 , 
1.7 x 10 -4 and 1.9 x 10 -4 s -1 , respectively. The 
code includes 80 rotational levels in the ground state 
( E < 943 cm -1 ) and in each of the vibrational states. 

Figure 5 shows how the fractional populations of 
the HDO rotational levels evolve from thermal to flu- 
orescence equilibrium. Calculations were performed 
with model parameters pertaining to comet C/1996 
B2 (Hyakutake) at it closest approach from the Earth 
(r^ = 1 AU, Qh 2 o = 2 x 10 29 s -1 ). Collisions with 
electrons were not considered in this calculation and 
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Fig. 5. — Evolution of the population of the rotational levels of HDO as a function of distance to nucleus. Calcu- 
lations are performed with Qh 2 o = 2 x 10 29 s -1 , r/, = 1 AU, v exp = 0.75 km s o c — 5 x 10 cm and 
= 70 K. Collisions with electrons are not considered. 


a c was taken to be 5 x 10 -14 cm 2 . Thermal equilib- 
rium breaks down at typically 500-1000 km from the 
nucleus for the low-lying rotational levels shown in 
the figure. This corresponds to the distance at which 
spontaneous emission rates begin to exceed collision 
rates. Radiative lifetimes are a few hundred seconds 
or less and the molecules relax rapidly to populate 
the lower states. At fluorescence equilibrium ( r h = 
1 AU), only the 0oo and l 0 i are significantly popu- 
lated, with 90% of the total population distributed 
among these two levels. This is because at r* = 1 
AU, the IR excitation rate that pumps the loi level 
is comparable to its spontaneous emission rate (Aij 
= 1.747 xlO -4 s’ 1 ). Between two fluorescence cy- 
cles, the HDO molecule has time to cascade down to 
the l 0 i rotational level, but not to fully relax down 
to the 0oo level. A fully relaxed population distribu- 


tion would be obtained at large heliocentric distances, 
while, well below 1 AU, a large number of levels would 
be significantly populated at fluorescence equilibrium. 

In Table 6, we give predicted brightness temper- 
atures for a selection of HDO millimeter and sub- 
millimeter transitions observable from ground-based 
telescopes. This may help for planning observations. 
We only consider two favorable cases for detection: 
1) a moderately bright comet at 1 AU from the Sun 
making a close approach to Earth (case of C/1996 
B2 (Hyakutake)); 2) an active comet at 1 AU from 
the Sun and from the Earth (case of C/1995 01 
(Hale-Bopp)). Because in both cases the observa- 
tions are sensitive to the intermediate region where 
neither LTE nor fluorescence equilibrium apply, we 
also give brightness temperatures calculated in these 
two limited assumptions. Calculations without elec- 
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tronic collisions are also given for comparison. 
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Table 1 

Log of the observations. 


Date 

[UT] 

Th 

[AU] 

A 

[AU] 

Int. time* 
[min] 

Frequency 13 

[GHz] 

Tsys 

[K] 

96 03 23.5 


0.118 

19 

464.925 

1500 

96 03 24.5 

1.06 

0.105 

31 

464.925 

1700 

96 03 26.5 


0.109 

49 

362.046 

1200 


integration time on source. 

^Central frequency of the main side band. 


Table 2 

Characteristics of the CSO telescope. 


Frequency 

[GHz] 


HPBW* 

["] 


Vmb h 

465 

362 


15. 

20. 


0.5 

0.7 

a Half Power Beam Width. 





b Main beam efficiency. 


Table 3 





Line intensities 



Molecule 

line 

Frequency 

[GHz] 

date 

[UT] 

/T B dt>‘ 
[K km s" 1 ] 

HDO 

loi~0oo 

464.92452 

96 03 23.5 
96 03 24.5 
96 03 23.5 - 24.5 

0.89 ± 0.16 
0.47 ± 0.15 
0.65 ± 0.10 

DCN 

5-4 

362.04648 

96 03 26.5 

< 0.18 b 


‘Main beam brightness temperature scale. The intensities and uncertainties are calculated on the high-resolution spectra over a 3 
km s' 1 velocity interval. 
b 3 <r upper limit. 
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Table 4 

Coma temperature measurements 


Species 

date 

Lines 

FOV 

[km] 

Trot * 

[K] 

Ref. b 

CO 

March 24-25 

346, 461, 691 GHz 

800-1600 

73 ± 6 

(1) 

CH3OH 

March 23-24 

462, 465 GHz 

1200 

69 ± 10 

(1) 

CH3OH 

March 24 

242 GHz multiplet 

1600 

75 ± 7 

(2) 

CH3OH 

March 30 

252 GHz multiplet 

4200 

65 ± 7 

(1) 

a Rotational temperature. 





b (l) From CSO observations, Lis et al. 

(1997); (2) from IRAM Plateau de Bure observations, Biver (1997). 




Table 5 





HDO, H 2 0 production rates AND (D/H)h 2 0 RATIO 



Species 

Date 

Q 

Source 

Reference 




[s- 1 ] 




HDO 

March 23-24 

(1.20 ± 0.28) x 10 26 

CSO 

this work 


h 2 o 

March 24.5 

(2.3 ± 0.3) x 10 29 

IRTF/H 2 o 

Dello Russo et al. (1997) 


March 25.1, 26.1 

(2.0, 1.4) x 10 29 

Nangay/OH a 

Gerard et at. (1998) 


March 23.1, 25.1 

(1.6, 1.9) x 10 29 

IUE/OH 

Festou et at. (1997) 


March 23.5, 24.5, 25.2 

(1. 6-2.5, 1. 7-2.7, 1.0-1. 7) x 10 29b 

Lowell Obs./OH a 

Schleicher (1997) 


March 23, 26 

(1.8, 1.9) x 10 29 

H Ly-a 

Bertaux et al. (1998) 


March 23, 24 

(3.1, 2.8) x 10 29 

Olf 1 /)) 

Hicks and Fink (1997) 

(D/H)hjO 

March 23-24 

(29 ± 10) x 10~ 5 





a Assuming (0h 2 O)/ < 3oh) = 1-1* 

b The range of production rates reflects calculations with different Haser parameters. The lower values (Schleicher 1997) use the Haser 
parameters of Cochran and Schleicher (1993). The larger ones are those computed with OH and OH-parent scalelengths at r h = 1 AU 
of 1.05 X 10 5 and 5.77 x 10 4 km, respectively, and an OH expansion velocity of 0.75 km s“ x . These latter parameters assume that 
A60H is collisionally thermalized and follows the H 2 0 outflow, as expected in the region sampled by the - 0.5 ' diameter aperture. 
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Table 6 

Line intensities of some HDO rotational lines? 


Line u [GHz] fT B dv [K km s~ 1 ] b 





d = 30 m c 



d = 10 m c 



Qh 2 o = 2 

X 10 s 9 s~ 

\ Tk 

= 1 AU, A = 0.1 

AU 





lio-lll 

80.578 

5.15 

(4.63, 3.90, 1.04) 

10 -2 

1.94 

(1.38, 1.28, 0.34) 

io- 

-2 

3i2“221 

225.897 

6.03 

(6.06, 5.98, 0.05) 

10~ 2 

1.74 

(1.72, 2.02, 0.02) 

10' 

-2 

2ii-2i2 

241.562 

1.09 

(1.08, 1.01, 0.07) 

10" 1 

3.75 

(3.53, 3.42, 0.24) 

10' 

-2 

220“3i3 

266.161 

4.14 

(4.07, 5.63, 0.02) 

10- 2 

8.30 

(7.80, 19.1, 0.08) 

10' 

-3 

loi“0oo 

464.925 

1.38 

(1.38, 1.32, 4.59) 


5.44 

(5.71, 4.49, 15.6) 

10’ 

-1 

3l2”3i3 

481.780 

1.05 

(1.06, 1.03, 0.01) 

10" 1 

3.44 

(3.45, 3.52, 0.03) 

io- 

-2 

2 o 2 ~lll 

490.597 

3.63 

(3.64, 3.46, 1.47) 


1.35 

(1.38, 1.16, 0.49) 



llO”loi 

509.292 

12.2 

(12.2, 11.6, 3.09) 


4.56 

(4.44, 3.97, 1.06) 



2ii-2q2 

599.927 

12.8 

(12.8, 12.4, 0.87) 


4.59 

(4.55, 4.24, 0.30) 



£ 

o 

II 

x 10 30 s~ 

1 , r h ■ 

= 1 AU, A = 1 AU 






80.578 

1.11 

(0.94, 0.81, 0.22) 

10- 1 

2.96 

(2.01, 2.02, 0.54) 

io- 

-2 

3l2~221 

225.897 

1.31 

(1.32, 1.39, 0.01) 

10- 1 

2.88 

(2.89, 4.15, 0.04) 

io- 

-2 

2ii~2i2 

241.562 

2.56 

(2.53, 2.36, 0.17) 

10- 1 

7.09 

(6.65, 7.11, 0.50) 

io- 

-2 

220“3i3 

266.161 

7.63 

(7.42, 13.3, 0.05) 

10- 2 

1.23 

(1.16, 4.03, 0.02) 

io- 

-2 

loi~0oo 

464.925 

3.57 

(3.61, 3.18, 11.1) 


1.47 

(1.57, 1.01, 3.51) 



3i2“3i3 

481.780 

2.51 

(2.53, 2.50, 0.02) 

10" 1 

6.95 

(7.01, 7.99, 0.07) 

io- 

-2 

2 o 2 “lll 

490.597 

9.26 

(9.36, 8.39, 3.57) 


3.29 

(3.41, 2.63, 1.12) 




509.292 

30.7 

(30.5, 28.2, 7.52) 


10.9 

(10.2, 9.05, 2.41) 



2ii~2o2 

599.927 

32.0 

(32.0, 30.2, 2.13) 


10.5 

(10.3, 9.78, 0.69) 




a Assuming a (D/H)h 2 o ratio of 25 x 10 6 (Qhdo/Qh 2 0) = 50 x 10 5 ). 

b Calculations are performed with the comprehensive model which includes col- 
lisions with electrons. The three values in parenthesis are those computed without 
electronic collisions, under LTE assumptions and at fluorescence equilibrium, re- 
spectively. The model parameters are Tkin = 70 K, v exp = 0.75 km s” 1 and a c = 
5 x 10“ 14 cm 2 . 
c Antenna diameter. 


23 


